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We study the emission of scalar fields into the bulk from a six- dimensional rotating black hole 
pierced by a 3-brane. We determine the angular eigenvalues in the presence of finite brane tension 
by using the continued fraction method. The radial equation is integrated numerically, giving 
the absorption probability (graybody factor) in a wider frequency range than in the preexisting 
literature. We then compute the power and angular momentum emission spectra for different values 
of the rotation parameter and brane tension, and compare their relative behavior in detail. As is 
expected from the earlier result for a nonrotating black hole, the finite brane tension suppresses 
the emission rates. As the rotation parameter increases, the power spectra are reduced at low 
frequencies due to the smaller Hawking temperature and are enhanced at high frequencies due to 
superradiance. The angular momentum spectra are enhanced over the whole frequency range as 
the rotation parameter increases. The spectra and the amounts of energy and angular momentum 
radiated away into the bulk are thus determined by the interplay of these effects. 



PACS numbers; 04.50.-h, 04.70.Dy 

I. INTRODUCTION 

Braneworld models with large extra dimensions 0,0) Q 
bring us an interesting possibility to address the hierar- 
chy problem by lowering the fundamental scale of gravity 
down to order of TeV. It has been argued in the con- 
text of TeV scale gravity that mini black holes might be 
created through high-energy particle collision at future 
colliders [J- Much effort has been directed towards a 
theoretical understanding of the black hole formation at 
TeV energies (e.g., Q). After their production, the black 
holes will decay via Hawking radiation |6j|. This process 
provides a window to probe high-energy physics, gravity 
at small distances, and properties of extra dimensions, 
which motivates recent extensive studies on this topic. A 
nonexhaustive samp ling of the literature can be found in 

Refs. p?, 's', ^, 1^, [il [H, Q E [M [13, [lip, S mi, 

|22, 23, 24, 25, 26, 231 . For a review see RefHi^Most of 
the related work to date has ignored the effect of brane 
tension, treating black holes as "isolated" ones (see, how- 
ever, Refs. [1^ 113, [U for the effects of self-gravity of 
branes). 

It is in general very difficult to obtain a black hole 
solution localized on a brane with finite tension be- 
cause tension curves the brane as well as the bulk (cf. 
Refs. [H, [U, m, [1^). However, codimension-2 branes 
exceptionally allow for a simple construction of local- 
ized black holes thanks to their special property; start- 
ing from the Myers-Perry solution [s^l one rescales the 
polar angle around a symmetry axis as ^ — > Bij) and 
then the brane tension is proportional to the deficit an- 
gle 27r(l — B). In this way both nonrotating [s^] and 
rotating [s^ black holes on codimension-2 branes have 



* Email: 
t Email: 
t Email: 



tsutomu" at" gravity, phys.waseda.ac.jpl 



nozawa" at" gravity.phys. waseda.ac.jp | 
takamizu" at" gravity, phys.waseda.ac.jpl 



been constructed. Following the work of [38|, Hawking 
evaporatiori [40| and the quasi-normal modes for bulk 
scalars [l!], and fermions [43|] have been investigated 
in the nonrotating background, showing that the finite 
brane tension modifies the standard result derived as- 
suming negligible tension. 

In this paper, we shall consider a six-dimensional ro- 
tating black hole pierced by a tense 3-brane and discuss 
the emission of massless scalar fields into the bulk. We 
intend to shed light on the spin-down phase in the life 
of a black hole, which is often neglected in the literature 
but could be of some significance. (In fact, a rotating 
black hole does not necessarily spin-down to zero, but 
evolves toward a nonzero angular momentum (4^. |45||.) 
Ignoring the brane tension, very recently Creek et al. 
studied the emission of scalars in the bulk in a higher- 
dimensional rotating black hole background \2^. They 
employed matching techniques to obtain an analytic solu- 
tion to the scalar field equation, which is a good approx- 
imation in the low-energy (wr^ ^1) ^^^d slow- rotation 
{a/rh ^ 1) regime, where lo is the energy of the emitted 
particle, rt is the black hole horizon radius, and a is the 
rotation parameter. In the present paper, with the help 
of numerical computations we are able to handle the in- 
termediate regime {ujrh ^ 1 and a/r^ ^ 1), and thus we 
not only include the effect of the finite tension but also 
extend the range of validity of [2^ . 

This paper is organized as follows. In the next section 
we give a quick review of the rotating black hole solu- 
tion on a codimension-2 brane. In Sec. Ill we present 
separated equations of motion for a massless scalar field 
and determine angular eigenvalues in the presence of the 
deficit angle. Then in Sec. IV the radial equation is solved 
numerically to give the power and angular momentum 
emission spectra. Finally we summarize our conclusions 
in Sec. V. Appendix contains the analytic calculation of 
the absorption probability, which complements the nu- 
merical results presented in the main text. 
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II. A ROTATING BLACK HOLE ON A 
CODIMENSION TWO BRANE 

We begin with a brief review of the rotating black hole 
solution on a codiniension-2 brane. (For further detail 
see Refs. [H, [11].) The solution shares some properties 
with the Myers-Perry black hole [33| . We are considering 
the models with five spatial dimensions, and so the rota- 
tion group is iS'0(5). The number of Casimirs (i.e., the 
number of mutually commuting elements of the group) is 
equal to rank[S'0(5)] = 2. Hence, we have two axes of ro- 
tation associated with two angular momenta. However, 
in the present article we will be focusing on the special 
but simple case of a single rotation parameter with the 
angular momentum pointing along the brane. This is in- 
deed an interesting case from the phenomenological point 
of view, because the black hole formed by the collision of 
two particles confined to the brane will have a single ro- 
tation parameter. The exact metric that describes such 
a rotating black hole is given by [39] 



_ ( 1 _ di^ - ^ sin2 edtdcf, + Cdr^ 



p^de^ + sin^ 9 



fia 
2^ 



■ sm 



■ cos^ e (dx^ + sin^ xdV'^) , 
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The coordinate ranges are < < 7r/2, < x < tt, and 
< < 27r. The parameter B is related to the brane 
tension a as 



B = 1 



27rM4 ' 



(3) 



where is the six-dimensional fundamental scale. We 
assume that < i? < 1. When B — 1, the above metric 
reduces to the usual Myers-Perry solution with a single 
rotation parameter in six dimensions [37| . When B ^ 1 
the solution is asymptotically conical. The parameters 
/t (> 0) and a denote the specific mass and angular mo- 
mentum, respectively, related to the ADM mass and an- 
gular momentum of the black hole as 



MBH = 2A4M4/iS, Jbh 



1 



(4) 



where Ad = 27r('^+i)/2/r[(d + l)/2] is the area of a unit 
d-sphere. Note here that the effect of the deficit angle 
B is separated out from the definition of the area. The 
black hole horizon radius follows from A{rh) = 0. For 
later purpose it is convenient to define the dimensionless 
measure of the angular momentum a* := a/rh- Since the 
sign flip a — s- —a simply changes the direction of rotation, 
in what follows we will assume a > without any loss of 
generality. 



Note that A(r;i) = has a root for arbitrary a. 
This should be contrasted with the four-dimensional Kerr 
black hole, which has an upper bound on the rotation 
parameter. This point will become clearer by consider- 
ing the black hole temperature. The generator of the 
horizon is parallel to the Killing vector = (9/9t)^ -I- 
^hid/dcj))'^, and the surface gravity of the black hole is 
given by ~ — (V'^^'^)(V^^i/)/2. Here flh is the angular 
velocity of the horizon: 



n 



h • = 



(5) 



Then, the Hawking temperature is given by 



H 



K 

2^ 



Anil + al)rh' 



(6) 



and one sees that the temperature never attains zero. 
This means the black hole is manifestly nondegenerate 
for any value of a. 

An interesting effect of the brane tension perhaps one 
is immediately aware of is the rescaling of the gravita- 
tional scale Af, M^^s = B^^^M^, as is observed in 
Eq. (j4|). This modifies the scenario of the black hole pro- 
duction and evaporation as follows [H, HI, 4fl] • Suppose 
that a black hole is formed through a particle collision 
and its mass A/bh is fixed by the accelerator. The up- 
per limit on the impact parameter b is given by 6 < 2rfi. 

Since Vh ^ /-^'^tis^ fixed A/bh the horizon radius 
is larger than the naive estimate, rn ~ B~^/^ , which is 
called the "lightning rod" effect [sl, |3^ . In the present 
situation producing black holes will become easier due 
to this effect. Since the black hole formed in this way 
has the angular momentum J = bM^Yi/'^i we have J < 
"/max := MBufh, and so the angular momentum can also 
be larger than naively expected, Jmax ~ Z?"^/"^. Now the 
time scale r of mass loss is estimated from dA^Bn/dt ~ 
(horizon area)xr^ ^ ^''^h^ ^ yielding r ~ B~^/^. Thus 
the lifetime of the black hole will become longer. Simi- 
larly, the time scale t' of angular momentum loss is esti- 
mated roughly from dJBn/dt ~ r/idA^Bn/di. This gives 
t' ~ Z?^^/"^, and hence the time scale r' will also become 
larger. 

The above mentioned modifications are caused by the 
rescaling of the gravitational scale. In the rest of the 
paper we separate this effect and discuss the impact of 
finite brane tension on Hawking emission with the horizon 
radius r^ fixed. 

Let us finally remark that in models with extra di- 
mensions the bulk will be compact. Throughout the pa- 
per we assume that the horizon size of the black hole 
is smaller than the typical compactification scale. With 
this, the above solution provides a nice description of a 
six-dimensional black hole localized on a 3-brane. 
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III. SCALAR WAVE EQUATION IN A 
ROTATING BLACK HOLE BACKGROUND 



We now turn to discuss Hawking emission of massless 
scalar fields. To this end we solve the equation of motion 
for a massless scalar field ^I' in the rotating black hole 
background ([1]). The governing equation is given by 



1 



-9 



=5,,(V^.g^'^9.*) =0, 



(7) 



where \/—g = Bp^r^ sin 9 cos^ 9 sin x- We can separate 
the above equation by assuming the ansatz 



^ = e'"^'e'""^R{r)S{9)T{x)E{i^), 
where m = 0,±1,±2, ... The radial equation is 



(8) 



dr 



dr 



1] + a Ll> 



where 

K :— {r^ + a^)Lj — am, A 
and the angular equations are given by 
1 d / . „ o „ d 



i? = 0, (9) 



2amu, (10) 



sm a cos^ 
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2 / 
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sin^ X 
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-i]S, 



(11) 
(12) 

(13) 



Here 77, v, and fi are separation constants. The above 
equations of motion are the same as those derived in 
Ref. [Ill except that in the last equation (fT^ there ap- 
pears "i3." We immediately see that 



n 
B 



with n = 0,±l,±2,. 



(14) 



As is emphasized in Ref. [40l|, the eigenvalues are cou- 
pled and therefore all of the eigenvalues will be modified 
when B ^ 1. Since the eigenvalues v and r] will be depen- 
dent on \n\, we can restrict to n > without any loss of 
generality. Negative values of n can be treated similarly. 

We proceed to determine the eigenvalue v. Performing 
the change of the variable and function as a; = cos x, 
T := (1 - a;2)-*'/2T, we have 

(i-^'):?^^-2(M + i)a;ir 



da; 



dx 



+ {iy - fi){iy + fi+l)T = 0. (15) 
The boundary conditions are given by 
(i^-/^)(i/ + Ai+l), 



dx 



±- 



2(m+1) 



-T at 



±1. 



(16) 
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FIG. 1: Eigenvalue rj for selected values of {j,£,m,n) and 
B. Solid and dashed lines refer to the series expansion (|28p 
truncated at 7th order, while points and crosses show the 
numerical results without relying on the small-aoj expansion. 



As in [43], we first consider the special case in which 
B = 1/N with N — 1,2,... and hence /i — nN is an 
integer. In this case it is easy to find v = fi, n + 1, n + 
2, ... Thus the eigenvalue v can be written in terms of an 

integer j as 



v = 3+K, 0„ := — g— 



(17) 



where j > \n\ and the absolute value signs are inserted 
for clarity. Although Eq. (fT7|) was derived for B^^ = 
l,2,...^his result can be generalized to arbitrary B E 
(0, 1] [42|. We have confirmed numerically that Eq. (fTT]) 
indeed holds for noninteger values of B~^. 

When B = I, V is independent of n (and as a result r/ is 
also independent of n). For the tensional case, however, 
v depends on n (and hence does rj) . From Eq. (|17p we see 
that the brane tension increases the eigenvalue i' relative 
to the tensionless case. 

To determine the eigenvalue 77, we exploit the contin- 
ued fraction method developed originally by Leaver [46[. 
Following 01, we write S* = (1 - u'^)\"'\/^u''S with 
u = cos 9. Then, S obeys 

d2 



[i^ + l)- — - (|m| + l)u 



du 



+ {a^u?u^ + 77 - v{v 3) - |7n|(|m| + 2iy + 3)] 5 = 0, 



(18) 



subject to the regularity boundary conditions 



—5 = at u = 0, (19) 

du 

d a^u? + T]~v{y ^-i) - \m\{\m\ -|- 2i^ + 3) 

— — O — j j o 

du \m\ + 1 

at u = 1. (20) 
We seek for a series solution in the form of 



S = ^ay^. 

p=0 



(21) 
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Substituting this to Eq. p8|) we obtain the three-term 
recursion relation 

aoai+/?nao = 0, (22) 
(XpCLp+i + f3pap + jpQp^i = 0, p=l,2,..., (23) 



where 



a 



p 



2(p+l)(2i^ + 2p + 3), (24) 
77- (;/+ |m| +2p)(//+ |m| +2p + 3), (25) 
{aujf. (26) 



The continued fraction equation for the eigenvalue is then 
given by 



In order for the series to converge in the limit auj 0, 
we require that it has a finite number of terms. Imposing 
f3q = for some integer q > and identifying 2q = 
i — {j + \m\), we obtain 



77= (f + 6„)(^ + 3 + 6„) 



(29) 



for acu = 0. Here we have a restriction for the integer £: 



1 



(£-j-|m|) = 0,1,2,... 



(30) 



00 



aoTi 



= 0. 



ai72 



(27) 



P2 



"273 
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The boundary conditions are automatically satisfied by 
the series solution (1^ . 

We expand the eigenvalue ry in powers of aui around 
acu — 0: 



To determine the coefficients in (|28p it is convenient to 
use the qth inversion of (P7)) : 



aq-2lq-l aq+llq+2 

Pq-1 - -3 Pg+1 



(31) 



•q-2 



Pq+2 



^ = ^VpiaujY 

p=0 



(28) With some manipulation we find ryi = 773 = and 



Vo 

f]2 

m 



(^ + 5„)(£ + 3 + 6„), 
2e{e + 3) + 3 + 2j + 2,f - 2m^ + A{£ + j + 2)b„ + Abl 
[2ii + bn) + l][2{i + bn) + 5] 



{£-j-\m\){£ + j-\m\ + l + 2bn) 



{e-j-\m\~2)(£ + j-\m\-l + 2bn) 



4[2(£ + 6„) + 1]2 
ji-j- \m\ + 2)i£ + j- \m\ + 3 + 2b„) 
4[2(£ + M + 5]2 



'?2 



4[2(^ + 6„)-l] 
(^ - j - |m| + 4)(j + ^ - |m| + 5 + 26„) 
4[2{£ + 6„) + 7] 



V2 



(32) 
(33) 

(34) 



One can confirm that Eqs. ((32|) ~ ((34)) correctly reproduce 
the known result in the case of 6„ = [H, S^l Although 
the expressions are too lengthy, higher order coefficients 
can be obtained easily. 

One can instead integrate Eq. (Ilip numerically to de- 
termine 77. The numerical computation has an advan- 
tage that it does not rely on the small-aw expansion. 
In Fig. [1] we compare the series expansion truncated at 
7th order with the numerical result that is free from any 
approximation. It can be seen from this that the an- 
alytic result can reproduce the numerical computation 



^ The second order coefficient in [47| contains a sign error, which 
has been corrected in |25||. 



remarkably well even for au ~ 6. Therefore, in the next 
section we will safely use the analytic approximation for 
77. We thus avoid numerical determination of the eigen- 
values and so the rest of the problem simply reduces to 
solving the ordinary differential equation for the radial 
mode function. 



We remark here that while the brane tension does not 
affect the n = mode, the eigenvalue 77 (for n ^ 0) 
increases as B decreases (and so the tension increases), as 
can be seen from Fig.[TJ This behavior has been observed 
for a nonrotating black hole in [40] , and now it turns out 
that the similar thing generally holds for a rotating one. 
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IV. POWER AND ANGULAR MOMENTUM 
EMISSION SPECTRA 

We are going to solve the radial equation ([9]) to com- 
pute the emission spectra. To do so we first specify 
the asymptotic form of the solution close to the hori- 
zon and far away from it. In terms of the new coordi- 
nate defined by dr* := [(r^ -I- a^)/A]dr and the function 



$ :— r{r'^ + a^Y^'^R, Eq. © can be rewritten as 



— $ + I b - mn(r)f - V{r)\ $ = 0. 
art ^ ' 



where ^(r) := a j (r^ ^ (P') and 
A 



(35) 



V{r) 



(r2 + a2)2 
1 



^2 r 



rir^ -|- a^^l"^ drj 



2U/2 



(36) 



It is easy to see that V(r) « as r — > r/j and r oo. 
Keeping this in mind the asymptotic solutions are found 
to be 



Th ■■ 



R^A. 



(oo)< 



A. 



(oo) e 

out ^2 



where 

k :~ LO ~ mQh^ 
and y is the tortoise-like coordinate defined by 



y := 



■ rh In 



A(r) 



(37) 
(38) 

(39) 
(40) 



We choose the boundary condition A^^j = 0, i.e., we 
impose that no outgoing modes exist near the horizon. 
The absorption probability is then given by 



\<Ajimn I 



A 



(oo) 
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(oo) 



(41) 
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FIG. 2: Absorption probabilities for {j,£,m.,n) — (1,1,0,1) 
5 = 1, 0.4, and various a*. 
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FIG. 3: Superradiance modes for a, — 1.2 and 0.6. The 
absorption probabilities are negative for oj* < mflhrh, where 
Qhr-h 0.49 for o* = 1.2 and Qhrh ^ 0.44 for a. = 0.6. 
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FIG. 4: Power emission spectra for different values of a, and 
B. 



Note the explicit dependence on the angular eigenvalue n 
of the absorption probability. This is due to the nonzero 
brane tension {B 7^ 1). 

Using the series expansion of the eigenvalue rj (trun- 
cated at 7th order) , we numerically integrate Eq. ^ and 
compute the absorption probability. We defer to Ap- 
pendix the analytic calculation of the absorption prob- 
ability in the limit lu^, Lurh <C 1 and a* <C 1. The 
analytic and numerical results are found to be in good 
agreement in the regime <C 1 and a* <C 1. 

As was remarked above, both and i' become larger 
as the brane tension increases. This results in the en- 
hancement of the "potential" (|36p . which will reduce the 
absorption probability. A typical example of the absorp- 
tion probability (as a function of w*) is plotted for var- 
ious and B in Fig. ^ showing that the absorption 
probability indeed decreases with decreasing B. We can 
see in Fig. [3] that there appear superradiant modes with 
< u) < mflh, for which the absorption probability is 
negative [isj . 

From the absorption probability we compute the en- 
ergy and angular momentum emission rates. They are 
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FIG. 5: Angular momentum spectra for different values of a, 
and B. 



given by the formulas 

1 w 2 

d^J 1 ^ m 2 

d^;: = ^ e^/T„_iiA-^™i ' (43) 

where Tf/ and k were already defined in Eqs. ([6|) and ([39]), 
respectively. We summed up to = 5 modes in calcu- 
lating the above quantities. Our numerical results are 
summarized in Figs. [Hand [S] 

We find that the finite brane tension reduces the power 
and angular momentum emission spectra. This is an- 
ticipated from the behavior of the absorption proba- 
bility stated above. The power emission rates are re- 
duced in the low frequency regime as the rotation pa- 
rameter a* increases. This is because the Hawking tem- 
perature becomes lower with increasing a*. However, 
the power emission rates are enhanced at high frequen- 
cies. This should be caused by superradiance. The an- 
gular momentum emission rates are enhanced over the 
whole frequency range with increasing a* because the 
effects of superradiance win out. (In the previous esti- 
mates [3, [il, significant superradiance is observed 
as the angular momentum increases.) One can see from 
Fig. [5] that a large portion of the contribution to the 
power emission spectrum is coming from the £ = m 
modes, which show the superradiant behavior (except 
for £ ~ m ~ 0). We confirmed that the same is true 
for the angular momentum spectrum. As a* increases 
the amplitudes of these modes are enhanced, lifting up 
the total spectra. One can also see that the oscillatory 
behavior of the spectra is due to these modes, with each 
peak corresponding to each mode. However, the above 
things do not mean that only the £ ^ m modes are im- 
portant in determining the behavior of the spectra; the 
eigenvalues for the £ = m modes are independent of the 
conical deficit and the other modes (|n| > 1) determine 
the decrease in the spectra with decreasing B. 
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FIG. 6: Contributions from £ = m modes. Red (solid) line 
shows the total power emission spectrum depicted also in 
Fig. |4j while green (solid) line represents the sum of £ = m 
modes, each of which is shown by blue (dashed) line. 

Since the brane induced metric is the same as the ten- 
sionless case, the Hawking spectrum of brane-localized 
fields does not depend on the tension (for fixed r/i).^ 
The emission of brane scalar fields has been studied in 
Refs. [H [U, m HH, [H for rotating black holes in vari- 
ous dimensions, and the brane-to-bulk ratio of the energy 
emission rates has been discussed in [25i] . implying that 
the dominant channels are the brane-localized modes. 
Since the finite brane tension further suppresses the bulk 
field contributions, it is likely that brane-localized scalar 
emission dominates bulk emission also in the present 
case. Although bulk modes will not be observed directly, 
bulk emission is still important because it indirectly de- 
termines the amount of energy and angular momentum 
left for brane-localized emission. 



V. CONCLUDING REMARKS 

The black hole production in TeV scale gravity offers 
us a possible window to explore the presence of extra di- 
mensions. In this paper we have studied Hawking emis- 
sion of scalar fields into the bulk from a rotating black 
hole localized on a codimension-2 brane. The exact so- 
lution we used is a rare example in which we can treat 
self-gravity of the brane consistently in the context of 
brane-localized black holes [HI. This simple model en- 
ables us to elucidate how the brane tension modifies the 
Hawking spectra relative to the tensionless case. 

Assuming the separable ansatz for the scalar field, we 
have determined the angular eigenvalue in the series ex- 
pansion form: rj = fjo + fi2{auj)'^ This analytic 



^ Although this statement is true in the present model, it is not 
clear whether or not the induced metric is independent of the 
brane tension in more reaUstic situations with compact extra 
dimensions. 
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approximation (truncated at 7th order) was turned out 
to be in excellent agreement with the numerical result 
in the regime auj < 6. Using the analytic form of ry, 
we then integrated the radial equation numerically and 
computed the power and angular momentum emission 
spectra. Our finding is that the finite brane tension sup- 
presses the power and angular momentum spectra. We 
also showed that the power emission rates are reduced at 
low frequencies and enhanced at high frequencies as the 
rotation parameter a* increases. The suppression at low 
frequencies is due to the smaller Hawking temperature 
and the enhancement at high frequencies is caused by su- 
perradiance. The angular momentum emission rates are 
enhanced over the whole frequency range as increases 
due to superradiance. The spectra and the amounts of 
energy and angular momentum radiated away into the 
bulk are thus determined by the interplay of these ef- 
fects. To conclude, the brane tension plays an important 
role in the evaporation process in the life of a mini black 
hole. 



It is possible to obtain an analytic but approximate 
solution to the radial equation, as has been done re- 
cently in [m and is replicated in Appendix. The range 
of validity of this analytic approximation is restricted to 
<C 1 and a* ^ 1. Therefore, even in the tensionless 
case (5 = 1) our result is new, in that we have extended 
the range of validity of (25| by invoking the numerical ap- 
proach. We however exploited the analytic expression for 
77 to simplify numerical calculations. The eigenvalues ry 
obtained here by the continued fraction method are quite 
accurate even for a<,uj^ ^ 6, which allows us to compute 
the spectra for tj* > 1 and a* > 1. In order to explore 
the full regime extending to auj > 0(10), we need to solve 
both the radial and angular equations numerically, which 
is left to further investigation. 



In this paper we considered only the emission of scalar 
fields. It would be interesting to study the bulk emission 
of higher spin fields and determine the brane-to-bulk ra- 
tio of the energy and angular momentum emission rates. 
Recently, it has been reported that for fermion fields the 
bulk emission dominates the brane-localized emission in a 
six- or higher dimensional Schwarzschild background [27| . 
Therefore, investigating the effects of the finite brane ten- 
sion [4^ and the black hole rotation on the fermion emis- 
sion would be of particular interest. Another open issue 
is to clarify the spin-down evolution of the rotating black 
hole and the effect of brane tension on it. This process 
can be studied along the line of [3, [1^ . We plan to re- 
turn to this issue in the near future. Finally, it is fair to 
say that all of the results in this paper have been derived 
ignoring the compactification mechanism, which may af- 
fect both brane and bulk emissions. This point is worth 
exploring, though constructing brane-localized black hole 
solutions in compact space will be quite difficult. 
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APPENDIX A: ANALYTIC APPROXIMATION 
METHOD FOR a;. < 1 AND a, <1 

In this appendix we present an approximation method 
to obtain an analytic expression for the absorption prob- 
ability. Our result here simply generalizes that of [2^ 
to include the deficit angle (see also [III). The proce- 
dure is as follows: first we obtain the asymptotic solu- 
tions in the far-field and near- horizon regions, and then 
match the two solutions in the intermediate zone. The 
approximation is valid in the low-energy (w* ^1) and 
slow- rotation (a* ^ 1) regime. We check that in this 
range the analytic result agrees with the numerical one 
displayed in the main text. 

Let us first focus on the near- horizon zone (r ~ r^). It 
is convenient to work with a new radial variable defined 

by 



for which we have 



dr 



(1 - i)rA{r) 



(Al) 



(A2) 



with Air) := 3 -I- a? /r^ . Now the horizon is located at 
^ = 0, while the asymptotic infinity corresponds to ^ = 1. 
Near the horizon, the radial equation ^ reduces to 



(1-70^ 



if? 



^(l-eM2 (1„^)^2 



i? = 0, (A3) 



where 



:= 3 + al, := (1 + al)uj^ — a,m, 

j:=l-4al/Al, A* A + i^(i^ + 1)0^. 

Performing the transformation i? = ^"(1 - ^)^i^(f ) with 

.K* 



(A4) 



we obtain a hypergeometric differential equation 

dF 



dbF = 0, (A6) 
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FIG. 7: Absorption probabilities computed analytically and 
numerically. Solid lines follow from the analytic approxima- 
tion, while dashed lines denote our numerical result. 



where a := a + /9 + 7 — 1, 6 := a + /3, and c = 1 + 2a. 
Thus, the near-horizon solution is given in terms of the 
hypergeometric function by 

R=A^Cil~O^Fi&,b,c;0 (A7) 

+ A+r^il - £.fF{a - c + 1, 6 - 5 + 1, 2 - S; 0, 

where A± are integration constants. In the hmit r ^ 
0) we have 

R ~ A-C + A+C = A^e-'''y + A+e'''y, (A8) 

where y is defined ear her in the main text. Since we are 
imposing the boundary condition such that no outgoing 
wave is present at the horizon, we set A^ = 0. Thus we 
arrive at 

i?NH = at r^rh. (A9) 

One can check that the convergence condition for the 
hypergeometric function, Re[c — a — 6] > 0, is indeed 
satisfied. 

Now we extend the solution (|A9[) to go beyond the 
near-horizon zone. Using the formula 

T{d)T{d-a-b) 



F{a,b,c;0 = 



-(1-0 



r(s-5)r(a-6) 

xF{a,b,a + b- c+ 1; 1 - 

c-a-imna + b-d) 



mm 

xF{c — a, c — 6, c — a — 5 + 1; 1 — 
and taking the hmit r — > cx) — > 1) we obtain 

i?NH ^ Air-^^ + yl2r-3(2-7-/3) ^ (aiq) 

where 



A, := A_(l + a2) 



/3^3/3 r(c)r(c-Q-b) 

T{d-a)T{d-b)'' 



A, := + „2)2-,-,^3(2-.-,) r(c)r(a + 6- c) ^ 

^ ^ T{a)T{b) 

(All) 

The expression (jA10[) should be matched to the far-field 
solution which will be derived below. 



Let us go on to the far-field solution. In the far-field 
zone (r ^ r/i) the radial equation ([9|) reduces to 



dr^ r dr V r 



The solution is given by 



Rff — B\ - 



B- 



Yritor) 



(A12) 



(A13) 



where Jr (Yr) is the Bessel function of the first (second) 
kind, T :— 1/77-1-9/4, and B12 are the integration con- 
stants. This solution is in turn to be extended to the 
near-horizon zone. Taking the limit ^ 1 we get 



Rpp 



p3/2p( 



Bi /ujr\'^ 

r(r + i) VT/ 



p3/2, 



2 / 



.(A14) 



The two solutions (|A10[) and (|A14[) appear to have 
different powers in r, but by taking the low-energy (cj* <g; 
1) and slow- rotation (a* <^ 1) limit we are able to match 
them. Neglecting terms of order ujI, al, and a*w*, we 
have A* ~ 3, 7 ~ 1, and (3 ~ —{£ + bn)/3, leading to 
-3/3 ~ f + 6„ and -3(2 - 7 - /3) ~ -(^ + &„ + 3) in 
Eq. (jAlOp . where we used 77 = {£ + bn){l + bn + 3) + 
Oialuil). As for Eq. jJTH, we have t -3/2 ~ £ + 6„ 
and -(r -I- 3/2) ~ -{i + bn + 3). Thus in this hmit we 
achieve exact matching. 

The absorption probability can be expressed in terms 
of the ratio of the coefficients Br :— i?i/i?2. One finds 



Br 



-(2f+26„+3) 



(A15) 



In the limit r — 
Rpp 



rr^(T)r(a)r(b)r(5-a-6) 
T{d-d.)r{£-b)r{d. + b-d)' 

-> 00 the solution (IA13P can be written as 
At 



If 1 /l(oo)l: 

J.2 ' ^out ^2 



(A16) 



where 



,(00) 



(Bi + iB2 



A 



(00) 



{Bi-iB2)- 



„i(r7r/2+7r/4) 

V27ra; 

-j(r;r/2+7r/4) 



V2^ 



(A17) 
(A18) 



The absorption probability is given by \AY 



l^out /^in I ■ 
BrB*r > (B,. - 



In the low energy limit (w* ^ 1) we have 
- -B;)/i > 1. Hence, 



2z(^-^ 



(A19) 



Substituting (|A15|) to (|A19|) results in a lengthy expres- 
sion, but when m is not large we may expand the equa- 
tion with respect to a (<j; 1) in the low-energy and slow- 
rotation regime. Then, at leading order in a = — iif^/A^, 
one finds 
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r2(2/3 + 7 - 2)r2(l - /3)(2 - 2/3 - 7) sin2[7r(2/3 + 7)] 



A, {£ + 6„ + 3/2)r2(£ + bn + 3/2)r2(/3 + 7 - 1) sin'[^(/3 + 7)] 

I 



+ ai) 



2U/3 



2^+2fc„+3 



(A20) 



In Fig. [7] we compare this analytic result with our nu- 
merical calculation employed in the main text. It can 
be seen that the numerical result indeed agrees with the 
analytic one in the low-energy and slow-rotation regime. 
Since all the terms except on the right hand side of 
Eq. (|A20[) are positive, the sign of \A\'^ is controlled by 
-ftT* = (l+a2)(jj» — a,m — rii{l+al){Lo — mflh) . Therefore, 
the absorption probability is negative for superradiant 
modes (0 < a; < mQh)- 

In the low energy regime the j— £ = 'm = n = mode 
will be dominant, for which the absorption probability is 
given by 



\A 



0000 1 



3(3 + a2)(2-7) 



(A21) 



This does not depend on the deficit angle. 

To obtain the expression for the absorption cross 
section we extract the ingoing s-wave from the plane 
wave 



IC 



(higher multipole moments), 
(A22) 



where the square root in the denominator is the normal- 
ization factor. In order to determine /C, one integrates 
both sides of Eq. (|A22p over the 4-sphere with a deficit 
angle, and then, looking at the far region wr 3> 1, ex- 
tracts only the ingoing modes. Thus, we arrive at 



(A23) 



leading to the low energy absorption cross section 



"87r2 



ri{ri + a')B 



(l + «^) 



(l + a2/3)(2-7) 



(A24) 



In the static limit we have tJo = (horizon area) , which re- 
produces the general result for the spherically symmetric 
case l50l. 
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